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Based on recent experimental evidences that the electronic charge degrees of freedom plays 
an essential role in the spin-Peierls-like phase transition of NaV2 05, we first make the mapping 
of low-energy electronic states of the d—p model for NaV2 05 to the quarter-filled t—J ladder 
with anisotropic parameter values between legs and rungs, and then show that this anisotropic 
t—J ladder is in the Mott insulating state, of which lowest-energy states can be modeled by the 
one-dimensional Heisenberg antiferromagnet with the effective exchange interaction Jefs whose 
value is consistent with experimental estimates. We furthermore examine the coupling between 
the ladders as the trellis lattice model and show that the nearest-neighbor Coulomb repulsion 
on the zigzag-chain bonds can lead to the instability in the charge degrees of freedom of the 
ladders. 
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Recently it has been reportecfl^ that the phase tran- 
sition at Tc = 34 K in NaV205 is not a conventional 
spin-Peierls transition but rather a charge ordering into a 
geometrical configuration of charges with the spin excita- 
tion gap. Althoughptkcrc^ai-jistill^lot of controversies in 
experimental data^Blfl'B§l3 SB Ef this work certainly 
casts doubt on the present understanding of the basic 
electronic states of this system; i.e., electronic charge de- 
grees of freedom cannot be eliminated in the low-energy 
physics of the system. More specifically, we might have 
to start from a quarter-filled ladder-type model, rather 
than from the one-dimensional (ID) Heisenberg model 
with only the dimerizing lattice degrees of freedom, to 
discuss this phase transition. 

In this paper we consider this unconventional spin- 
Peierls system. We first analyze the low-energy states 
of the d—p model by using a V4O4 cluster and by mak- 
ing one-to-one correspondence of the states to analyt- 
ically solved states of the 4-site t—J cluster we derive 
the anisotropic t—J ladder Hamiltonian which we ar- 
gue can hopefully be a starting model for the low-energy 
physics of this system. We will show that the map- 
ping works quite well and enables us to evaluate the 
parameter values for the Hamiltonian. A Lanczos exact- 
diagonalization technique is then applied to the finite- 
size clusters of this ladder model to examine its low- 
energy states. We will thereby show that our anisotropic 
t—J ladder has a finite charge gap in the realistic param- 
eter region and the lowest-energy levels can be modeled 
by the ID Heisenberg antiferromagnet with the effective 
exchange interaction consistent with experimental esti- 
mates. We will moreover argue that because the model 
still has the charge degrees of freedom in its low-energy 
sector it can be a starting Hamiltonian for describing the 
unconventional spin-Peierls system NaV2 05 if appropri- 
ate lattice degrees of freedom and/or some presently un- 
known factors are taken into account; we will suggest a 



possible scenario that a fairly small inter-ladder repulsive 
interaction can induce instability for the charge dispro- 
portionation in the ladders. 

High-energy physics of the system may be contained 
in a d—p model cluster V4O4 (a pair of the rungs) which 
is illustrated in Fig. 1. The d—p Hamiltonian reads 
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where d\^ (pia) creates (annihilates) a spin-cr electron on 
the fia:j,-orbital at site i (neighboring p-orbital at site I), 
and nf^ and n^^ are the electron number operator. We 
take into account (i) the hopping parameters tpd between 
the dxy-orhital on V ion and the p^j-orbital (t^^) ^^'^ Py' 
orbital (i^^) on the neighboring O ions where the phase 
factor (±)' has to be taken appropriately, (ii) energy- 
level difference Ap^ = ef — between the da;j,-orbital 
and the Pa;-orbital (A^^) and p^-orbital (A^^), and (in) 
Coulomb repulsion Ud on V ions. The anisotropy of the 
parameter values will turn out to be important. We note 
that this cluster is the smallest one in that it takes into 
account the anisotropic legs and rungs simultaneously; 
i.e., it contains the effective repulsive interaction on the 
rungs which can lead to the charge gap, the effective 
exchange interaction between the rungs which can lead 
to the ID Heisenberg antiferromagnet, and the exchange 
interaction between the neighboring V ions on the leg 
of the ladder which can lead to the spin gap when the 
system is dimerized at T < Tc- The coupling between 
ladders is neglected for the moment (see below for this 
effect). Values of the d—p model parameters are taken 
from Ref.y): t^^ and t^^ take the value of 1.22 and 1.03 
eV, respectively, but are |assumed to depend on the bond 
length d as tpd « ^^7/2 ^) ^ g 5 ^y, A^^ = 4 eV, 
and Ud — 4: eV. There are 10 electrons in the cluster. 
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Fig. 1. Orbital structure of (a) tiie V2O3 ladder in NaV2 05 and 
(b) a corresponding t—J ladder. The small clusters used for the 
mapping are indicated by dotted lines. 



Energies of several low- lying states of the V4O4 clus- 
ter are evaluated by a Lanczos diagonalization method 
and are compared with the analytically solved low-energy 
states of the 4-site t—J cluster with two electrons. The 
t—J Hamiltonian reads 

Ht-j = - 5] ty(cLc,,+H.c.) 

+ ^J,,(s,.S,-in,n,) (0.2) 

where (ij) represents the nearest-neighbor bonds along 
the legs and rungs with the hopping and exchange pa- 
rameters tij and Jij taking t± and J± for the rungs and 



t\\ and J|| for the legs, c, 



„t 



(1 — rii^a-) is the con- 



strained electron-creation operator at site i and spin a 
(=Ti J.) J Si is the spin-i operator, and rii ~ ni| + riii is 
the electron-number operator. We confirm that the sev- 
eral low-lying states of both d—p and t—J clusters have 
the same quantum numbers and characters and thus we 
find that the mapping is straightforward. Here we de- 
termine the values of the 4 ladder-parameters by fitting 
the energies of the lowest 5 states for the d—p and t—J 
clusters (the extra one is for the origin of energy). The 
obtained values of the ladder parameters for the realistic 
bond lengths are: t±_ = 0.298, t\i = 0.140, J± = 0.049, 
and J|| = 0.029 in units of eV (which we will use in 
the following discussions unless otherwise indicated). In 
Fig. 2, we show the bond-length dependence of these pa- 
rameters. We note that with increasing dy , the value of 
J|| decreases as expected, but the value J± also decreases 
considerably. Also noted is that with increasing d± the 
value of J|| increases. These are from the higher-order 
contributions to the exchange interactions. The hopping 
parameters t± and <|| on the other hand vary more or 
less as expected. 

The single-particle spectra from the quarter-filled 
ground states of the V4O4 cluster and corresponding t—J 
cluster are shown in Fig. 3. We find a very good agree- 



Fig. 2. Bond-length dependence of the estimated parameter val- 
ues for the t—J ladder. In (a) and (b) the V-O bond length 
on the rungs d± is kept constant at d± = 1.825 A, while in (c) 
and (d) the V-O bond length on the legs dy is kept constant at 
dn = 1.916 A. 



ment, demonstrating that the mapping is satisfactory 
also for the one-electron added and one-electron removed 
sectors of the cluster Hilbert space. 

Now let us examine the electronic states of the ob- 
tained anisotropic t—J ladder at quarter filling. Here we 
examine the charge gap Ac of the model, which is nat 
yet well known although there were some discussions. 
The charge gap is defined by 
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[Eo{Nt + 1,N^)-Eo{Nt,N^)] 
-[EoiN^,N^)-EoiN^-l,N_^)]\ (0.3) 



where Eo{Ni, Ni) is the ground-state energy of the sys- 
tem with 7V| up-spin and down-spin electrons. This 
is evaluated for 8, 12, and 16-site clusters with either 
periodic or antipcriodic boundary condition to minimize 
the shell effect on the values of Ac. The results are plot- 
ted in Fig. 4 as a function of 1/L {L is the length of 
the ladder) to see the value at the infinite system size. 
We find that at least for t_L/t|| > 1.3 — 1.5 the gap opens 
up although the situation is still unclear for t± ~ t|| in 
the present calculations of small finite-size systems. We 
also calculate the equal-time charge correlation function 
{niUj) and find that the probability of finding two elec- 
trons on the same rung is very small for t±/t\\ > 1.3, in- 
dicating that one electron is localized on each rung. It 
seems thus quite reasonable to say that the anisotropic 
ladder at quarter filling is in the insulating state for the 
realistic parameter values. 

The mechanism of this charge localization is common 
in a class of quarter-filled organic systems where the 
Mott insulating phase is realized due to dimerization: 
with increasing dimerization strength a dimer turns into 
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Fig. 3. Single-particle spectra for the d—p cluster compared with 
those of the t—J cluster at quarter filling. The spectra near 
the chemical potential (set ai u) = 0) are enlarged in the lower 
panel. For the d—p model, the total spectra summing up the 
contributions from the V d^y-orhital and O Pa;- and p^-orbitals 
are shown although around the chemical potential the spectral 
weight is almost completely from the V d^^y-orbital. 



Fig. 4. Charge gap Ac/t|| of the finite-size clusters for the 
quarter-filled t—J ladder as a function of the inverse of the clus- 
ter size 1/L. Anisotropy in the hopping parameters f^/ty is 
varied with keeping Jx/J|| = 2 and Jx/tx = 0.2 constant. 
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an effective single site of a half-filled antiba: 
with Hubbard-like repulsive interaction.ta-llj^li 
mechanism operates also in the anisotropic ladder at 
quarter filling if wc can regard the rung as a dimer, and 
actually in some parameter region, we do have this in- 
sulating phase as we have seen above. This Mott insu- 



lating phase can be modeled as a ID Heisenberg antifer- 
romagnet of which the effective exchange interaction is 
estimated to be Joff = 85 meV from the singlet-triplet 
splitting in the 4-site t—J cluster. This-value is con- 
sistent with the value obtained in Rcf.EiP and should 
also be consistent with values determined from the ob- 
served temperature dependence of the uniform suscepti- 
bility x{T^) if ^ small inter-chain exchange ^cnupling and 
some other effects are taken into account. El Note that 
the value of Jy is much smaller than the value of J^s 
(see Fig. 3); the picture that the electrons order on one 
of the legs of the ladder already above Tc should pre- 
dict a very different peak position in x{T). Also noted 
is that the bond-length dependence of Jeff is very strong 
and anisotropic (see Fig. 3), which may be confirmed 
by some high-pressure experiment. The resultant pic- 
ture for the high-tcmpcrature phase of J^aV205 is thus 
equivalent to that of Horsch and MackElP although the 
procedure of derivation of the ID Heisenberg model is 
quite different. 

In the real materials, we have to take into account 
the coupling between the ladders, i.e., the system may 
be the one modeled as a two-dimensional (2D) trellis 
lattice (ladders connected with zigzag chain bonds as in 
Fig. 1 (b)). We examine the electronic states, i.e., the 
equal-time charge correlation function and charge gap 
for the 16-site clusters of this lattice (or coupled two 
4x2 ladders), and find that the basic electronic state 
that one electron is localized on each rung is maintained 
against the inter-ladder hopping integral of up to txy — 
1.2<||. The reported values of t^y are fairly small, so 
that the electronic statCLnLthe real materials may be in 
this localization regime. ErliiP We note that the situation 
is similar to the case in the 2D dimerized lattice system 
of K-(BEDT-TTF)2X where the Mott insulatinftJitate is 
realized over a wide range of parameter values .EZP 

We now want to consider the low-temperature phase 
of NaV205, i.e., the mechanism of charge ordering re- 
ported in Ref.Er. This is in principle possible because 
our starting Hamiltonian, the t—J ladder at quarter fill- 
ing, still has the charge degrees of freedom in the low- 
energy—mgions, unlike in the coupled Heisenberg chain 
in Ref . We however immediately notice that it seems 
unlikely that the t—J ladder with the small J/t values 
has inherent charge instability in itself even though it is 
anisotropic; the only possibility seems to be the phase 
separation realized for very large J|| values where the 
electrons are localized on one of the legs to form a Heisen- 
berg chain. 

The simplest interaction that can lead to the charge 
instability is the intersite Coulomb repulsion. In 
(TMTTF)2X, e.g., the nearest-neighbor Coulomb repul- 
sion is discussed to_be the origin of the observed charge 
disproportionation.E^ Here we introduce the repulsive 
interaction V between ladders, i.e., on the zigzag chain 
bonds running along the crystallographic 6-axis, which 
have the shortest bond-length among V-V bonds. We 
use the two 4x2-ladder clusters for modeling the trel- 
lis lattice and calculate the charge correlation function 
{uirij). As shown in Fig. 5, the tendency to the charge 
disproportionation where the electrons are localized on 
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one side of the legs is actually found for not too large val- 
ues of V (only 10-20% of Ua)- Once the ID Heisenberg 
chain on one of the legs of the ladder is thus realized (by 
some phase transition in the thermodynamic limit) the 
system may be subject to the lattice dimerization which 
leads to the opening of the spin excitation gap. 



Fig. 5. Charge correlation function (nirij) for the t—J treUis lat- 
tice as a function of the density repulsion V between the ladders. 
The site indeces defined in Fig. 1 (b) are used with i = 1 . The pa- 
rameter values for the zigzag chain bond arc assumed as txy = tii 



and Jx 



: J|, 



An optimistic speculation for the phase transition of 
NaV205 is a realization of this charge ordering. In this 
picture, the critical temperature of the charge ordering 
needs not be the same as the spin-Peierls transition tem- 
perature of the spin chain of one of the legs of the ladder; 
if the former temperature is higher than the latter, we 
should observe two transitions independently. However 
the opposite seems to be the case in real materials: when 
the charge ordering occurs, the spin chain of one of the 
legs of the ladder is already below its spin-Peierls tran- 
sition temperature, so that the chain dimerizes simulta- 
neously at Tc into the dimerization geometry potentially 
prepared at a higher temperature. A significant devia- 
tion from the BCS-type temperature dependence of the 
gap functiond' might be attributed to this.EJ' Also inter- 
esting is a report from mean-field calculations of the ID 
extended Hubbard model for (TMTTF)2X where it is 
shown that the first-order phase transition of the charge 
ordering crosses over to a second-order-like behavior 
when the ratio oLjthe values of two hopping parameters 
deviates from l,!!^ which is suggestive to the observed 
second-order-like phase transition of NaV2O5.0' How- 
ever, to decide if the mechanism suggested above really 
applies, we of course have to know much about the elec- 
tronic states of this system from further experimental 
studies. 

In summary, we have studied the electronic states 



of NaV2 05 by making the mapping of low-energy elec- 
tronic states of the d—p model for this compound to the 
quarter-filled t—J ladder with anisotropic parameter val- 
ues between legs and rungs. We have shown that this 
anisotropic t—J ladder is in the Mott insulating state, of 
which lowest-energy states can be modeled by the one- 
dimensional Heisenberg chain with the effective exchange 
interaction Jeff whose value is consistent with experi- 
mental estimates. We have also studied the coupling 
between the ladders by using the trellis lattice model 
and have shown that the intersite Coulomb repulsion on 
the zigzag-chain bonds can leads to the instability in the 
charge degrees of freedom of the ladders. A possible 
mechanism of the spin-Peierls-like phase transition of 
NaV205 has thereby been suggested. 
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